Current fluctuations of an interacting quantum dot 
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We calculate the counting statistics of electron transfer through an open quantum dot with charg- 
ing interaction. A dot that is connected to leads by two single-channel quantum point contacts in 
an in-plane magnetic field is described by a Luttinger liquid with impurity at the Toulouse point. 
We find that the fiuctuations of the current through this conductor exhibit distinctive interaction ef- 
fects. Fluctuations saturate at high voltages, while the mean current increases linearly with the bias 
voltage. All cumulants higher than the second one reach at large bias a temperature independent 
limit. 
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The statistical distribution of current fluctuations in 
nanoscale conductors, the so-called "counting statistics", 
has received considerable attention over the past decade 
- theoretically [H 0] as well as experimentally U ■ The 
theory of this distribution is well developed for non- 
interacting conductors. Effects of electron-electron inter- 
actions, however, have only been studied in a few limiting 
cases so far. Interacting conductors in the tunneling limit 
of weak transmission have been shown to exhibit Poisso- 
nian current fluctuations 0. Also for weakly interact- 
ing conductors at arbitrary transmission the statistics of 
charge transfer has been found to be qualitatively identi- 
cal to that of a non-interacting system - it is multinomial 
i i. In a conductor that has arbitrary transmission 
combined with strong interactions one expects more pro- 
found changes of that statistics. This is indicated by a 
classical analysis that has already revealed qualitative 
changes of the distribution of current fluctuations in that 
regime. In this Letter, we confirm this expectation by 
a nonperturbative quantum mechanical analysis of non- 
Gaussian current fluctuations in a strongly interacting 
nanostructure at arbitrary transmission. 

An interacting nanostructure is at low energy scales 
described by the model of a non-interacting conductor in 
an electromagnetic environment, that is a quantum con- 
ductor with an effective series resistor 8]. Recently, Safi 
and Saleur have shown, that for a single-channel con- 
ductor this model maps onto a Luttinger liquid (LL) with 
impurity. The statistics of fluctuations for this model has 
been obtained at zero temperature in the context of tun- 
neling in fractional quantum Hall samples |lO| . From a 
perturbative analysis one expects that the current / in 
systems describe d by a LL scales with a power of the ap- 
plied voltage V [llj. In the limit of a strong impurity, 
one has / oc l/^/f"^, where g is the interaction parame- 
ter of the LL. The current Ib that is backscattered by a 
weak barrier obeys Ib csc V^^~^. Evidently then, g — 1/2 
is a special point in parameter space: The backscattered 
current becomes independent of the applied voltage. By 
virtue of the mapping of Ref. , this unusual behavior 
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FIG. 1: QD connected to two leads by single-channel point 
contacts QPC 1 and QPC 2. QPC 2 introduces backscatter- 
ing. The current / of charge transferred from lead 1 to lead 
2 is measured. 



should be observable in a quantum conductor with series 
resistance I/Gq, where Go = e^/h is the conductance 
quantum. An open quantum point contact (QPC) is a 
natural realization of such a series resistor. 

Motivated by this, we study transport through the in- 
teracting quantum dot (QD) shown in Fig. ^ It is con- 
nected to two leads via single-channel QPCs. One of 
them (QPC 1) is open such as to act as the desired series 
resistor to a scatterer in QPC 2. An in-plane magnetic 
field allows only electrons in one particular spin state to 
pass the QPCs. We address the incoherent case when 
inelastic relaxation inside the QD is faster than electron 
escape. The charging energy of the QD is assumed to 
be much larger than the temperature T and the applied 
voltage. In our analysis, we employ techniques that have 
been developed in the context of the Coulomb blockade 
1^ 13. 14, 15.J and charge pumping in almost open QDs 
[l6|. After bosonization, the QD in Fig. ^ is at low ener- 
gies, indeed, described by a LL at g = 1/2 (the Toulouse 
point). This allows us to study its transport behavior 
nonperturbatively by refermionization. We find as ex- 
pected that Ib, the current backscattered by QPC 2, 
saturates with increasing voltage. It is, however, hid- 
den in the total current I = Iq — Ib under a large back- 
ground current Iq = gGoV of charge carriers that are not 
backscattered. Analyzing also the fluctuations of the cur- 
rent through the structure, we find that the same physics 
leads to a saturation with increasing V of non-Gaussian 
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current fluctuations as well. No similar background as Iq 
is present for these fluctuations, since all electrons that 
are not scattered produce only Gaussian thermal noise. 
This, together with recent experimental advances in the 
detection of non-Gaussian current fluctuations brings 
the observation of this counterintuitive saturation effect 
within experimental reach. Remarkably, all cumulants 
higher than the second one are entirely temperature inde- 
pendent in the limit V oo. This is in marked contrast 
with the behavior of non-interacting conductors. 

We assume the QD large enough for electrons to re- 
lax inelastically before they escape the dot. In this case, 
the QD can be modeled by one one-dimensional electron 
mode for each QPC subject to a charging interac- 

tion inside the QD (at a: > 0, cf. Fig.^). At the energies 
e of interest, that are smaller than the Zeeman energy 
of electrons in the applied magnetic fleld, only electrons 
of one particular spin projection contribute to transport 
through the QPCs. For typical conductors the electronic 
spectrum at the relevant energies can be linearized and 
the Hamiltonian is diagonalized by bosonization. The 
QD is then described by two bosonic fields 9j and (f)j for 
every lead j that describe electron density fluctuations 
|l4l Il7j. Moreover, for energies below the charging en- 
ergy Ec of the QD, e Ec, the "charge" mode 6i +62 is 
pinned and the system is described by an effective Hamil- 
tonian iJeff for the "transport" modes dg — {di ^ 62) /2 
and (ps — {(f>i — </>2)/2. Os{0) is the operator of the charge 
Q transported through the QD, Q = e6's(0)/7r. We have 

{h = kB = 1) 
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A is the backscattering strength of QPC 2. 

To quantify the current / through the QD and its fluc- 
tuations, we compute a generating function Z that gen- 
erates moments of the charge Qr = /J" dtl{t) transferred 
during time t. 
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where rCk is the fcth cumulant of the distribution. Writ- 
ing Z as a Keldysh path integral [l^ and integrating out 
(j)s as well as all modes 9s{x) at x 7^ 0, we arrive at 



:,+Ss-i J (it2A(cos2e+-cos2e") 



(3) 



with the free action 
duj 



So = 



Lo [ANe+*e- - {2N + i)(|6'+p + \e-\^)] 



for fields 9^ corresponding to the mode 9s{0) on the for- 
ward (0+) and the backward {d~) part of the Keldysh 



time contour. N = [exp(aj/T) — 1] ^ is the Bose-Einstein 
distribution. The source term 



Ss = i I dt 
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couples V and ^ to charge eO/n and current —edt9/Tr, 
respectively We eliminate Ss by the change of in- 

tegration variables 0^{t) 0^{t) — eVt/A± e^/8, where 
V — V — 2i^T. The resulting Z can equivalently be 
obtained by integrating out all modes 4>{x) and 9{x) at 
a; ^ in a path integral corresponding to the Hamilto- 
nian expression 



if^dtH-(t) -if^dtH+(t) 
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with time-dependent Hamiltonians 



(5) 
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T± orders operators along the Keldysh contour and 

„2 



2o(C) = exp 



(6) 



(7) 



is the generating function in the absence of backscat- 
tering. We now follow the standard procedure 0| to 
refermionize H^. For this we define new fields 



ix)TH-x)]+9ix)±9{-x) (8) 



and a chiral Fermion b — a^^^^ exp(i(^_|_) with a short 
distance cutoff a. then have a representation as non- 
interacting Hamiltonians for b. Following Matveev 0| 
we introduce a Majorana fermion 0? -f and define new 
fermion operators c by the relation b = (d + d^)c. This 
brings into a quadratic form in fermion operators. 



eVt/2 



}■ (9) 



Inserting Fermion coherent states, we rewrite Eq. jSjl in 
refermionized form as a path integral |20|. We again 
integrate out all modes c{x 7^ 0) and are left with an 
integral over vector fields c = (c+(0), c~(0)) and d = 
(d+,d-). 



d*G^'d 



m = 2o(0 j -DcVd cxp{-* j ^ 
+ cG-\ + V^A {[d + d*)e'^^^'I^T'c + /i.e.)] }, 



(10) 
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where is the third Pauh matrix. Gd{oj) = 
—i{T±d{ijj)S{uj)) = —t'- jio is the Green function of 
d and that of c corresponding to the Hamiltonian 
© at A = 0. The time dependence of the scatter- 
ing term has been removed by a gauge transformation 
c exp(jeW/2) c that shifts the frequency of Gc- Since 
Gc is an electron Green function Gc{x,x') evaluated at 
coinciding spatial coordinates x — x' — 0, it is linearly 
related to the semiclassical Keldysh Green function of c 



1 - 2/ 2/ 
2(1-/) 2/-1 



(11) 



Gc = Gs/^ivF- Here, f{uj) = [cxp(cj - eV/2)/T + 1]-^ 
is the Fermi distribution function after the gauge trans- 
formation. 

The action in Eq. H10(l is diagonal in frequency and Z 
consequently factorizes into contributions from different 
frequencies. The Gaussian integrals result in 



lnZ(0 =lnZo(0+r 



duj 
2^ 



ln<^ 1 



rf2 



i)/+(i-/-) + (e-'^«-i)r(i 



(12) 



with /+ = /, /-(c^) = 1 - /(-c^), and Tb = j^vp. 

We remark that one arrives at an equivalent result by 
applying the quasiparticle formalism developed by Fend- 
ley, Ludwig, and Saleur [iJl to a LL at 5 = 1/2. 

From Eq. (O we find for the first three cumulants 
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where ip is the digamma function. At zero temperature, 
Eqs. (|13|l - H15|l obey the known relations between the 
higher order cumulants and Ci derived in Ref. ^3 • 

We first analyze the limit of a large voltage eV ^ T, Tb 
(while eV ^ Ec). In this limit, the backscatterer is weak. 
Charge transfer is then best understood by singling out 
two contributions: First, transfer of electrons that are 
not backscattered. They are responsible for Iq. Second, 
the contribution due to backscattered electrons, generat- 
ing Ib which converges to a limiting value — {e/A)TB 




FIG. 2: Cj, is plotted in units of eGoTb versus eV/Ts for 
different temperatures. For large voltages, C3 saturates at 
the temperature independent value — —{■Ke/8)GoTB ■ The 
inset shows that, at zero temperature, C3 changes sign at an 
effective transmission FeH « 0.26. 



[ll| . In the total current I ~ Ci, Ib is, however, hid- 
den under the contribution Iq — GoV/2 of unscattered 
electrons that increases linearly with V. Similarly, the 
large voltage hmit = Gq{ttTb/4: + T) of the second 
cumulant, Eq. (|14|l . has these two contributions: the fluc- 
tuations due to the backscatterer as well as the Johnson- 
Nyquist noise GqT of electrons that are not backscat- 
tered. Non-Gaussian current fluctuations, in contrast, 
allow to specifically probe the backscatterer and its large 
voltage behavior. This is because unscattered electrons 
produce purely Gaussian fluctuations. Accordingly we 
find that all higher order cumulants saturate at large 
voltages and, remarkably, their limiting values are tem- 
perature independent (see Fig. El for C3). The generating 
function in this limit takes the form 



lnZ°°(0-r 



Go 



2 



(»• 



ei/2 



(16) 

The temperature independence of higher order cumulants 
in the large voltage limit is in stark contrast to the behav- 
ior of a non- interacting conductor P|. It is thus a clear 
signature of interactions in the QD considered here. It 
can be understood by noting that a large voltage bias 
shifts the energies at which electron occupation num- 
bers are thermally smeared far away from the equilibrium 
Fermi level. Electrons with these high energies are, how- 
ever, effectively not backscattered and therefore produce 
purely Gaussian noise. Thermal fluctuations do thus not 
contribute to higher order cumulants. With present day 
experimental techniques Q this anomalous behavior can 
be observed in C3, the third cumulant of current fluctua- 
tions, Eq. (|15|) . It has the temperature independent large 
voltage limit Cf = —{en/8)GoTB, as shown in Fig.|21 In 
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the opposite limit of small voltages, we find at low tem- 
peratures T <^ eV <^Tb (moderately low such that the 
inelastic processes that justify our model are still opera- 
tive) the statistics 



lnZ(0 = 




Due to the scaling of the impurity strength with energy, 
the high and the low voltage limits, Eqs. H16|l and (|17|l . 
correspond to the weak and the strong backscattering 
limit, respectively. They are Poissonian, in accordance 
with Refs. 0,^3- 1°'^ voltages, as manifest in Eq. 
(|17|l . charge is transferred in units of the elementary 
charge. The current scales as / cx F'^, as expected from 
perturbative calculations 0|. The high voltage statis- 
tics, Eq. l(T^ . suggests, that charge is transmitted in 
packets carrying half the elementary charge. This is a 
direct consequence of the charging interaction of the QD 
that after every backscattering of an electron induces a 
positive electric potential on the QD. It attracts electrons 
from the leads to compensate for the electron that is miss- 
ing on the QD. In response, an electron flows onto the 
QD from either lead 1 or from lead 2. This either cancels 
or completes the transfer of an electron through the QD 
that was initiated by the backscattering event. Both pro- 
cesses are equally likely in the weak backscattering limit 
of an almost open contact QPC 2. Therefore only ev- 
ery other backscattering event transfers charge through 
the QD. Equivalently one can say that every such event 
transfers only e/2, as indicated by Eq. l(TB|) . 

In the shot noise limit of zero temperature, the third 
cumulant for non-interacting electrons vanishes at trans- 
mission r = 1/2. This is because in this case every elec- 
tron is transmitted or reflected with the same probability 
1 /2 independently of all other electrons. The distribution 
of transferred charge is consequently symmetric around 
its mean and its skewness C3 vanishes. This intuition re- 
mains correct for weakly interacting conductors j^. One 
expects it to be invalidated, however, by strong inter- 
actions that correlate the transfers of different electrons 
with each other. The inset of Fig. |21 shows that this is, 
indeed, the case for the QD we consider. We define an 
effective single-electron transmission Fcff — I/{GoV — I) 
for QPC 2 in series with QPC 1 by Ohm's law. C3 van- 
ishes then for T^s ~ 0.26, in clear contrast with a non- 
interacting structure. 

In conclusion, we have studied current fluctuations in 
a strongly interacting quantum conductor at arbitrary 
transmission. While it had been found in Refs. 0, 
that interactions do not qualitatively change the statis- 
tics of current fluctuations in perturbative situations, 
our nonperturbative solution does display features that 
are qualitatively different from those of non-interacting 
structures. This makes the measurement of current cor- 
relations a promising tool to probe interactions in the QD 
we considered and most probably in many other strongly 



interacting conductors. More specifically, we find that 
the fluctuations of the current saturate at high voltages, 
while the current itself increases linearly with the applied 
voltage. Moreover, all cumulants higher than the sec- 
ond one are temperature independent in the high voltage 
limit. We have discussed in detail how the voltage and 
temperature dependence of the third cumulant display 
these and other qualitative interaction effects. Experi- 
mental techniques for its measurement are available 01 ■ 
Our predictions arc thus experimentally testable. 
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